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Abstract 

This paper presents an adaptive tracking control method for a class of nonlinearly parameterized MIMO dynamic systems with 
time-varying delay and unknown nonlinear dead-zone inputs. A new high dimensional integral Lyapunov-Krasovskii functional 
is introduced for the adaptive controller to guarantee global stability of the considered systems and also ensure convergence of 
the tracking errors to the origin. The proposed method provides an alternative to existing methods used for MIMO time-delay 
systems with dead-zone nonlinearities. 
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1 Introduction 

In the past decades, control of multiple-input-multiple- 
output (MIMO) practical systems has attracted 
a great deal o f attention such as flyin g robots 
( Ge. Ren. Tee, and Lee Id2009l) ). b iped robots (Ge. Li. Yangl 
(20111)) and underwater vehicle (ICui. Ge. How. Chop I 
( 20101) '). To characterize certain non-sensitivity for 
small control inputs of these MIMO systems, dead-zone 
nonlinearities have to be considered otherwise it can 
severely limit system performances, even leading to 
instability. To handle dead-zone nonlinearities, many 
control appro aches have b een d evelop ed in the litera¬ 
ture s uch as dZhang. Ge I (120071) ) and dChen. Liu. Lin I 
(120131) ). On the other hand, time-delay usually can be 
encountered in thes e MIMO systems, dealt wi th by 
many authors, e.g., (ICui. Ge. How. Chop I (l2010P ) and 
references therein. Time-delay, like the dead-zone non- 
linearities, can also degrade the system performances 
or lead to instability if ignored during the course of 
controller designs. Given the effects of the dead-zone 
and time-delay, this paper considers a class of nonlin- 
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early parameterized MIMO dynamic systems where 
time-varying delay and unknown nonlinear dead-zone 
are simultaneously taken into accounts. 

To the considered problem in this p ap er, the most rel¬ 
evant pa pers are dZhang. Gel (120071)). dChen. Liu. Linl 
( 20131)). (IZhoul (I200S)). (IShvu. Liu. Hsu I (12005)) an d 


( Hua. Wang. Guanl (|2008l) ). In dZhang. Gel (I2007D ). 

an adaptive neural control was proposed for a class of 
uncertain MIMO nonlinear state time-varying delay 
systems with unknown nonlinear dead-zones and gain 
signs, however, it only guarante es semiglobal stability o f 
the c losed-l o op sy stem. Both dChen. Liu. Lin I d2013l) ) 
and (|Zhou I (|2008l) ) used the backstepping techniques 
to construct the controllers. However, the backstepping 
procedure is computationally time-consuming due to the 
computation of many virt ual controllers for the M IMO 
systems. The authors of dShvu. Liu. Hsu I d2005l) ) pro¬ 
posed a decentralized controller for large-scale systems 
wit h time-delay and dead- zone nonlinearities. However, 
in (|Shvu. Liu. Hsu I d2005l) ). the time-delay is constant 
and the parameters of th e dead-zone are known. In 
dHua, Wang, Guanl d2008l) ) the authors constructed a 
novel Lyapunov function and then designed a smooth 
adaptive state feedback controller for the time-delay 
system with dead-zone input. However, linear dead-zone 
was considered and the tracking error can converge only 
within an adjustable region. Thus one may wonder if it 
is possible to propose a new control method completely 
different from those existing ones to overcome these 
disadvantages mentioned above? This paper provides 
an affirmative answer by introducing a new high di¬ 
mensional integral functional as a Lyapunov-Krasovskii 
































































function of the closed-loop systems. 

From the motivatio n above and following our previous 
work (|Ge. Li I (l2014f) l. this paper introduces a new high 
dimensional integral Lyapunov-Krasovskii functional 
for a class of nonlinearly parameterized MIMO sys¬ 
tems with time-varying delay in states and unknown 
nonlinear dead-zones to achieve its tracking control. 

Compared to the existing results, the main contribu¬ 
tions of this paper are: i) By proposing a Lyapunov- 
based adaptive control structure, neither cancelation of 
the coupling matrix during linearizing the system nor 
conventional backstepping techniques is needed; ii) By 
introducing a new high-dimensional integral Lyapunov 
function in the control design, the process of controller 
design is simplified, i.e., it is unnecessary to calculate 
the inverse of the unknown control gain matrix; iii) By 
the construction of the Lyapunov-Krasovskii functional, 
the unknown time-varying delay in the upper bounding 
function of the Lyapunov functional derivative can be 
easily eliminated; iv) The developed control strategy is 
applied to a 2-DOF robotic manipulator system and the 
comparative simulation studies demonstrate the superi¬ 
ority of the proposed method. 

2 Problem Statement and Assumptions 

Consider the following uncertain MIMO nonlinear time- 
delay system with dead-zone nonlinearities 

X, — , 

x„ = B -1 (x) [F(x(t - r(t))) + D (t) + u], (1) 

y = xi, 


eigenvalues. Therefore, for the positive definite case, we 
have the following inequalities 

A m m(B)||a;|| 2 < x T Bx < A rora (B)||x|| 2 , Vx £ IT, 

where X m in{B) and A maa: (B) denote the minimum and 
maxim u m eigenvalues of B. respectively (Li. Ge. Wang] 
12008\ )). (Li. Dina. Gao. Duan. Su\ (2013 )). 

Assumption 2.1 (Zhang, Ge I ( 2001 )) The dead-zone 
outputs Hi,... ,u m are not available and the dead-zone 
parameters bi r and bn are unknown constants, but their 
signs are known, i.e., bi r > 0 and bn < 0, i = 1 , 2, ..., m. 
The growth of the ith dead-zone’s left and right functions, 
gn(vi) and gi r (vi), are smooth, and there exist unknown 
positive constants ku o, km, ki r o and ki r \ such that 

0 < kao < g a (vi) < km, Vn, £ (-00,6^], (2) 

0 £ ki r 0 £ 9ir(.Ui) — ^irl, 'dVi £ [6j r ,TOo), (3) 

where /3io < mm{kuo, ki r o} is a known positive con¬ 
stant, and g u (vi) = dgn(z)/dz\ z = Vi and g ir {vi) = 
dg ir (z)/dz\ z = Vi . 

Assumption 2.2 (Zhang, Ge I (200A) ) The unknown 
state time-varying delays Ti(t) satisfy fi(t) < f max < 1, 
i = 1, • • • ,m, with the known constants fmax- 

We know that there exist (IZhang, Ge I (120071 )1. fu{vi) £ 
(—oo, bi r ) and fi r {vi) £ (bu, + 00 ) such that 

9it(vi) = gu(vi) - gu(bu ) = gu{.fii{vi)){vi - bn), 

for fii{vi) £ (vi,bu) or ( b a ,Vi ), (4) 

Qir{Vi ) — 9ir(Vi) 9ir{bir) — 9ir((lir{Vi)){Vi bi r ), 

for £,ir(Vi) £ (Vi,b ir ) or (b ir ,Vi). (5) 


where x = [xf, x|\ ..., x^] T £ is the state vec¬ 
tor, Xj = [xa,x i2 ,..., x im ] T £ K m and = Xj(i = 

1,2,..., n); the nonlinear function F(x(t — r(t))) £ R m ; 
B (x) £ K mxm are unknown continuous bounded func¬ 
tion matrix; and the nonlinear function vector D(f) = 
[d\,d 2 - ■ ■ ■, dm\ T denotes the external disturbance, u = 
[u\,U 2 , ■ ■ ■ ,« m ] T £ is the output of the dead-zone 
control input and satisfies that if Vi > bi r , Ui = gi r {vi)', 
if bn < Vi < b ir , Ui = 0; if Vi < b iU m = gn{vi), 
where Vi £ R is the input to the ith dead-zone, and bu 
and bir are the unknown parameters of the ith dead- 
zone. In the paper, we consider the dimension of y and 
Xj(i = 1,2,..., n) are equal, therefore, y = Xi holds. 


Define vectors = [(p ir (t), (pu(t)] T and Ki(t) = 

L Sir (6r (Vi (t))), g'u Hu (' Vi(t)))} T with 

1 if Vi(t) > bu, j 1 if Vi(t) < b^, 

<Pil{t) = < (6) 

0 if Viit) < bu, { 0 if Vi(t) > b ir . 

Based on Assumption 2.1, the dead-zone control input 
can be rewritten as follows 

u i = K I (t)®i(t)Vi + £i(Vi), (7) 

where 


The control objective is to find a control input 
V = [vi,V 2 , ■ ■ ■, v m ] T such that the output of the system 
tracks the desired trajectory yu £ R m , while all the 
signals of the closed-loop system are globally bounded. 

Remark 2.1 The introduction of B(x) is for the phys¬ 
ical meaning of mechanical system, which can represent 
to the inertia matrix for the system. 

Remark 2.2 The matrix B(x) £ R mxm is known to be 
either uniformly symmetric positive definite or uniformly 
symmetric negative definite for all x £ R”, and have m 


Qir (fi r (i-i ) ) bir, If Ui £ bi r , 

-[o'iMuivi)) 

+9ir(£ir{Vi))]Vi, if bu < Vi < b ir , 
k -g'u{f,ii(vi))bu, if Vi < bu, 


( 8 ) 


and |£i(i’i)| < p*, where p* is an unknown positive con¬ 
stant with p* = {k^ 1 + kui) m&x{bir, —bu}. Therefore, 
we have 

u = K t <$>V + 3, 


(9) 



























where K = diag [Kfit)], $ = diag[$j(t)] and 3 = 

Let Brf(x) £ R mxm be a diagonal matrix with diagonal 
elements bdu 0(i = 1,2,..., m), then, there exists an 
unknown matrix Ab such that B(x) = B^(x) + Ab is 
satisfied. We can rewrite (1) as 

B(x)x n = F(x(t — r(i)))+D(t) + u. (10) 

Substituting (9) and B(x) into (10), we can obtain 
B d (x)x„ = T(x(t - r(t))) + K t $V + A P , (11) 

where 

T{yi{t - r(i))) = (I - A b B' 1 (x))F( x(i - r(i))) 

= [/i,/ 2 ,...,/ m ] T (12) 

and 

A P = -A B B- 1 A: T d>V + (I - A b B- 1 )(3 + D) (13) 
are column vectors. 

Assumption 2.3 I (7 c. Hana. Zhana\ fl99fi) ) Func¬ 
tions baa and di are continuous unknown, bdu and di re¬ 
spectively satisfy b d u = anddi < W di $ di (t) 

where Wbu £ R ; is unknown bounded constant param¬ 
eter vectors, <hpjj(x) £ M 1 is the known continuous 
smooth bounded regressor vector, W d i £M. 1 is a vector of 
unknown bounded constant parameters, and $*(1) £ R z 
is a vector of smooth bounded nonlinear function x. 

Assumption 2.4 KZhana. Ge I \200S ID The unknown 
continuous functions f\, fi,..., f m satisfy the inequality 

m 

\fi\ < Afc( x fc(i - T fc (t))) (14) 

k= 1 

with Qik{xk{t)){i = 1,2, ...,m) being known positive 
continuous functions. 

Lemma 2.1 iGe. Lee. Harris I 1(1991 j) ) LetH(s) denote 
an (n x m)-dimensional exponentially stable transfer 
function, r be the input and e = H(s)r be the output. 
Then r £ L ™ f) L™ indicates that e, e £ Ltf f) e 
continuous, and e —>• 0 as t —> oo. Moreover, if r ^ 0 as 
t —> oo, then e —» 0. 

3 Control Design and Analysis 

Define the filtered tracking error Si (ISlotinc. Li I (ll993ll ) 

Si = e[ ^ + Xne[ ■*+... + Xi^n-ia, (15) 

where e» = - y d i, e-" _1) (i = 1, 2,..., m) is the (n - 

l)th derivative of e*, A,i,..., A^ n _i are positive con¬ 
stants and are appropriately chosen coefficient vectors 
such that a —> 0 as s* —» 0 (i.e. r m_1 + Xnr m ~ 2 + ... + 
Xi >n -i is Hurwitz). 


From (15), we have 

S = B ( ^ 1 (x) [.F(x(i - r(i))) + K T <S>V + A P ] + v, (16) 

where S = [si,..., s»,... , s m ] T and v = [iq, ..., zq,..., 
V m ] T with 

v% = -Vdi + Aiie^ ^ + ■ • ■ + Aj in _i. (17) 

We now construct a new high dimensional Lyapunov- 
Krasovskii functional(see Eq. (42)). The first part of the 
Lyapunov-Krasovskii functional is chosen as 


V x = S T BaS, 


where 

1 

' 1 

= / 'd'B a d'd = diag 

f ’dB aii (5ti)d'd 

J 

0 

J 

.0 


with B q = B d a = diag \b d uOin]mx.m and matrix a £ 

For easy analysis, we choose an = ... = a mm . By 
exchanging x n i in x with dsi + (i(i = 1,2, ...,m), 
we define x, = [xf, x(f,..., x^, x nl , x n2 ,. .., iUst + 

£ R” m where Q = y* _1) - & with 

= Ajie^” 2 ' ) + ... + Xi yU -iei. i9 is a scalar and inde¬ 
pendent of x,;. We can choose suitable B^(x) and a, 
such that bductu > 0. 

Remark 3.1 We propose B^, which is a diagonal ma¬ 
trix, on the basis of weighted control Lyapunov function 
(WCLF) defined in (Ge, Hang, Zhang (1999)). The di¬ 
agonal element of B^ is defined as i9B ai i(xi)di9. By 
introducing (19), we construct (18) as the first part of the 
novel high dimensional Lyapunov-Krasovskii functional. 
Equation (18) is positive definite in the filtered error S 
and grows unbounded as ||S|| —> oo. 

Because B# in (18) depends on time t , the time deriva¬ 
tive of Vi includes the differentiation of matrix B d with 
regard to time t. T o facilitate com putation of its deriva¬ 
tive, according to ( Gentle (2007?)), we introduce a ma¬ 
trix operator for derivative operation of matrix-value 
function with respect to time t , i.e., for a time-dependent 
matrix A £ R mxn and a vector b(t) £ M ( , a matrix op¬ 
erator A4g(A, b) £ R mxn is defined with the entry of its 

*th row and jth column being Aidij(A , b ) = b with 

i = 1,2,... ,m and j = 1,2,... ,n. 

Differentiating (18) with respect to t gives 

Vi = 2 S t B^S + S T A4 a (B. 5 , S)S 

+S r Af a (B fl ,x)S + S T Af 3 (B,,C)S, (20) 

where < = [Ci, C 2 , • • •, tm] T , M d (B^,S) £ R mxm , 
A4 a (B tf ,x) G R mxm and A4g(B d , () £ K mxm are given 






















below 


Since matrices B d , a and B d a are symmetric,we have 


Ala(Btf, S) = diag 


Xa(B^,x) = diag 


C) = diag 


LO 


OSi 


nm otd 

* s 

3=1, 




( 21 ) 

( 22 ) 

* = 1,..., m(23) 


Let a = $Sj(i = 1,2,...,m), we can obtain 

9B aii <!B a 2 i do clB a 2 j 

9sj <9cr 9sj 9 ct 

•9B a 2j 9B aii 9(j 9Bo,j2 

9$ 9ct 9$ 9ct S *’ 


(24) 

(25) 


Noting that 9 is a scalar and independent of Q, and the 
fact Q = —Vi(i = 1,2,..., m), we have 

(26) 

0 0 


Motivated by (21)-(26), the following equations can be 
obtained 


S T M a (B rf ,S)S = S T [d 2 B Q ] — 2 / dB a d9 S 


= S t B q S-2S t BoS, 

S J M a (B tf ,()S = S T ( - f ] u 


(27) 


= -S 2 B a v + S 2 / B a i/d#. 


(28) 


By using the two equations above, we can rewrite (20) as 


Vi = S t B q S - S T B a v 

l 

+S T A4a(B^,x)S + J B a vd& 


(29) 


Using (16), we have 

Vi = S t B q B" 1 (x) [T(x(t - r(i))) + i^ T $U + A P ] 


+S 1 


A4 a (B^,x)S+ / B al'd.'d 


(30) 


B a B d 1 (x) = B d (x)aB d 1 (x) = a. 

Then, equation (30) could be rewritten as 
Vi = S T a [j-(x(t - r(t))) + K t $V + A P ] 


+S J 


A4 a (B^,x)S + J B a ud/d 
o 


(31) 


(32) 


From Assumption 2.3, we can rewrite (32) as 

Vi = S T a[J r (x(t - r(i))) + W T A> + K T <f>V + A P ],(33) 

where 2 = [x T , x T , S T , v T ] T , W = diag[W Bii ], 4> B = 
diag[<h Bil ] and 


1 1 
^'(x) = J t?A4g(<h B ,x)Sd$ + J $ B i/di?. 


(34) 


Then, given Assumption 2.4 and Lemma 2.1 in (iGe. Li I 
(1201411 ). it is easy to rewrite (33) as follows 


Vi = S T a W T A>(z) + K t $V + A : 

+S T aJ r (x(t — r(t))) 

< S T a W T ^(z) + K t $V + ||S T a||(7i + 72 

m m 

+yl|S T a|| 2 + - T fe (t))). (35) 

j=l k =1 


Remark 3.2 According to Remarks 2.1 and 2.2, we can 
obtain that the bounded matrix B(x) is either uniformly 
symmetric positive definite or uniformly symmetric neg¬ 
ative definite for all x £ K". Then the inverse B _ 1 (x) 
and unknown matrix A p should also be bounded. Accord¬ 
ing to Assumptions 2.1 and Assumptions 2.3, S, K and 
the disturbance vector D (t) are bounded. Therefore, the 
vectors (/— AbB _1 (x))E, (/-AbB _1 (x))D(1) andma- 
trix AbB _1 (x) are bounded. As a result, there exist two 
positive parameters 71 and 72 such that 

||(/-A B B- 1 (x))(2 + D(t))||<7i, 

|| - AbB-^xJ^VH <72||U||. 

Simila r to literatures iGe. Li\ \20lA )) and KXu. Ioannou\ 
US;;, 71 and 72 are regarded as robust parameters in 
controller (36). In practice, we can constantly adjust the 
two parameters until the proposed controller can stabilize 
the closed loop system. 

It is easy to check that 4t(^) is well-defined even if S 
approaches zero. We design an adaptive control 

V = A ssn ui + u 2 , (36) 

ur = -/V 1 ((*i + y<*)S M + Wfi*|.|(*) + T h ) , 


(37) 









































(46) 



P = 


71+72 ||ui || 



if S ± 0, 
if S = 0, 


(38) 

(39) 


where A sgn = diag[sg , n(sj)]; (*)|.| denotes matrix or vec¬ 
tor that its every element is the absolute value of (*)’s 
corresponding element; pis positive when 72 < ||/9jj" 1 || 1 ; 
W is the estimate of W; K\ is a positive diagonal ma¬ 
trix; T will be defined later; and /3q = diag[/3io, • ■ •, /3 m o] 
where /3io has been defined in Assumption 2.1. 

Since > /3o > 0, when S ^ 0, we can obtain 

S T aK T <f>V < -S T a(ATi + ^a )S - Sj r | a(IT | ' z j’I'|.| (*) 

+T M )-||aS||p 

< — S T a(Ki + — a)S — S T a(W T (z) 

+T) - ||aS||p, (40) 


From the definition of V a , we have 

m 

i=1 

The time derivative of V\j. (t) is 

-j m 

v 1 > ) w = 2 (,- f 

^\ ‘max) i L 
-&jk( x k(t - T k (t)))(l - Tfc(t)) 


Thus, the time derivative of V 2 is 


V 2 =V 1 +V a + ^Vj / .(t) 
. 7=1 

< S T a(—KiS — T) 

1 


2(1 7~max ) 


EE^^W)- 


3=1 k=1 


(47) 


(48) 


where we use the following facts 

—Sj r |Q;W’| 7 j4'j.|(2) < — 

-S^aT|.| < -S T aT. 

Since \\V\\ < ||ui|| + ||u 2 ||, we have ||S T a||( 7 i+ 7 2 ||l / ||)- 
||aS||p < 0. 

Substituting (40) into (35) and noting that ||y|| < 
l|ui|| + ||u 2 || give 


where T = [Ti,..., X*,..., T m ] T with 


T,; = 


2(1 T max )Si(Xii 


^Qik^kit)). 


k= 1 


(49) 


Noting that if T is utilized to construct the con¬ 
trol law, controller singularity may occur, since 
(2(1 - f max )s t a ti )~ 1 J2k=i e%k( x k(t)) is not well-defined 
at Si = 0. Therefore, define T,; as follows: 


Vi < -S T aiFiS - S T aW T 'S>(z) - S T aX 

-.mm 

+^EE Qjk( x k(t - T k (t))). 

3=1 fe= 1 


(41) 


T 

2 


2(i-fl x ) SM ELI if s i ± 

0 , if Si = 0. 


(50) 


Then, when Si ^ 0, 


Consider the following Lyapunov function candidate as 

m 

V 2 =V 1 +V a + ]Tv c/ .(f), (42) 

1=1 

m 1 

vE^W- ( 43 ) 

2=1 Z 

where Wi = Wi — Wi. 


V 2 < -S T aATiS < 0. (51) 

To this end, our main result can be summarized as: 

Theorem 3.1 For the closed-loop system (1) and (36), 
under Assumptions 2.1-2-4, f or bounded initial condi¬ 
tions, the tracking error e\ converges to zero, and the 
overall closed-loop control system is globally stable in the 
sense that all of the signals in the closed-loop system are 
globally bounded. 


The adaption law is designed as 

^ = n^i{z)a iiSi , (44) 

where 11, > 0 (i = 1,2, ...,m) is a diagonal constant 
matrix to be designed. 


Remark 3.3 From Lemma 2.1, it is obtained that the 
equation S = 0 defines a time-varying hyperplan in R n on 
which the tracking error e\ converges to zero asymptoti¬ 
cally. On the basis of this conclusion, we can obtain that 
e\ converges to the origin if S converges to zero asymp¬ 
totically. 


V jj j (t) is introduced to overcome unknown time-delays 
Ti(t), T 2 (t),..., T m (t) and defined as 


1 m r 

= on _ - -vE] / P 2 jk( x k(T))dT. (45) 

Tmax) J 




Remark 3.4 In practice, to prevent chattering phenom¬ 
ena, sgn function in (38) should be replaced by sat func¬ 
tion, i.e., 


S/e s , if ||S|| < e s , 
S/1|S||, if ||S|| > e s , 









where e s is a small constant. 


cr W 2 = 0.1, a v \ = a v2 = 0.01, = T w2 = diag{2.5}, 

r„i = r„2 = diag{15} and f max = 0.6. 


4 Simulation Results 

To validate the proposed method, we consider the fol¬ 
lowing 2-DOF robotic manipulator system 


) xi = x 2 , 

x 2 = B d ~\x) [.F(x(i - r(t))) + D(£) + K T $V] , 
y = xi, 

where B ( j _1 (x)J r (x(£ — r(t))) = B _1 (x)F(x(£ — 

r(£))), B d (x) = diag[6 d n,6 d22 ], D(£) = [d 1 ,d 2 ] T , 
F(x(£ - r(£V)) = [oi ,a 2 ] r , B(x) = [bu,b 12 -,b 21 ,b 22 ], 
xi = [qi,q 2 ] T , x 2 = [qi,q 2 ] r , b dll = 2m 2 hl 2 (cosq 2 + 2), 
bd .22 = m 2 l 2 2 ( 1 + 0.5sin q{), ai = —m 2 hl 2 qi(t - 

T 2 )q 2 (t-T 2 )sinq 2 (t-T 1 )-m 2 l 1 l 2 q 2 (t-T 2 )(q 1 (t-r 2 ) + 
q 2 (,t-T 2 ))smq 2 {t-Ti), a 2 = m 2 l 1 l 2 ql(t - t 2 ) sin q 2 (t - 
n), &n = (mi + m 2 )l\ 2 + m 2 l 2 2 + 2m 2 lil 2 cosq 2 , 
b\ 2 = b 2 i = m 2 l 2 2 + m 2 hl 2 cosq 2 , b 22 = m 2 l 2 2 , 
d\ = (mi + m 2 )hg cos q 2 + m 2 l 2 g cos(qi + q 2 ) and 
d 2 = m 2 l 2 cos(gi + q 2 ). 

We choose robot parameters as mi = 6 kg, m 2 = 3 
kg, l\ = 0.8 m, l 2 = 0.4 m and g = 9.81 m/s 2 for nu¬ 
merical simulation. We consider the desired trajectory 
yd = [sin£, cost] T and set the initial conditions q(0) = 
[0.3,0.5] T and g(0) = [0,0] T . We choose Wi(0) = -0.1 
and W 2 (0) = 1 as the initial value of adaption law. The 
design parameters of the above controller are An = 10 , 
A 2 i = 5, 7i = 7.2, y, = 0.2, fU = 0.02, 0 2 = 0.1, 
K\ = diag[10,5], /3o = diag[0.3,1.1], a = diag[14,14] 
and e s = 0.2. The parameters of the dead-zone are given 
as Qiriyi) — fair (.Vi bi r ) and gu{Vi) — foil (.Vi bn'j 
with the parameters of the dead-zones ku = 0.5, k± r = 
1.5, k 2 i = 1.5, k 2r = 2.5, bn = -0.5, b lr = 0.5, b 2L = 
—2.5 and b 2r = 2. The time-varying delays Ti(t) = 
0.2(1.1 + sint), t 2 (£) = 1 — 0.5cos t and f max = 0.6. 

The tracking errors between the joint positions and their 
references are shown in Figs. 1-2. The time histories of 
the adaptive parameters are shown in Fig. 3. These three 
figures show good transient performances the proposed 
method achieves. 

To show the advantages of the proposed method, we 
choose the conventional P D contro l a nd the adaptive 
neural control proposed in flZhang. Gc I (j2007|)) for com¬ 
parisons under the same time-varying delays and un¬ 
known dead-zones. As a traditional control method, the 
PD control can be written as V = —K p e — K d e. In the 
comparison simulation study, K p and K d are respec¬ 
tively set to K p = diag[45, 35] and K d = diag jl O, 1 0]. 
For th e adaptive neural control proposed in (IZhang. Ge I 
(|2007f) bwe also use two 3-layer neural networks contain¬ 
ing 10 hidden no des to approximat e the unknown func¬ 
tions as done in (IZhang. Ge I (I2007ID . The controller pa¬ 
rameters are chosen as 710 = 720 = 3.5, c Sl = c S2 = 0.1, 
An = 7.5, A 2 i = 5, ku = 3, k 2 1 = 3.2, A73 = k 2 3 = 0.02, 
V 1 = 0.15, 72 = 0.2, pi = p 2 = 2, o\ — fj 2 = 0.01, a wl = 


The comparative simulation results are shown in Figs. 1 
2 and Figs. 4-5. From Figs. 1-2, we can see that the PD 
method cannot make the tracking errors converge. From 
Figs. 4-5, the control sig nals of the adaptive neural con¬ 
trol proposed in (IZhang. Ge I (l2007f) j can cause chatter¬ 
ing phenomenon, which in practice can degrade system 
performances. From the simulation results, our method 
can have better performances. Moreover, we construct 
the controller mathematically by using the adaptive 
technique to deal with the uncertainty of the considered 
system, i nstead of us i ng ne ural n etworks approxim a- 
tion as in (IZhang. Ge I (12007IP and (IZhang Ge 1(120081) ). 
Furthermore, the method of (IZhang. Ge I (l2008h 4 lies in 
the backstepping technique which needs to construct n 
controllers in the n steps while our method only needs 
one step in the sense of backstepping. 


In order to investigate the control performances of the 
proposed method under different controller parameters, 
we also choose different parameters in the simulation. 
Specifically, we select three pairs of different values of An 
and A 2 i, i.e., case 1: An = 10, A 2 i = 5, case 2: An = 5, 
A 2 i = 2.5, case 3: An = 2.5, A 2 i = 1.25, to observe how 
these two parameters affect the control performances. 
Figs. 6-7 and Figs. 8-9 show the tracking error trajec¬ 
tories and controller output trajectories under different 
controller parameters, respectively. From these figures, 
we can observe that the greater the values of An and A 2 i 
are, the faster the convergence rate of tracking errors is, 
but accordingly the larger the control signals are at the 
beginning of t = 0. 


5 Conclusions 

In this paper, an adaptive control method is proposed 
for a class of uncertain MIMO nonlinear time-varying 
delay systems with unknown nonlinear dead-zones. The 
design is based on the use of a new high dimensional 
integral functional as a Lyapunov-Krasovskii function of 
the closed-loop systems, with the advantages of global 
stability and convergence of the tracking errors to origin. 
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